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ABSTRACT 

We study the stellar populations, star formation histories and star formation prop- 
erties for a sample of blue compact dwarf galaxies (BCDs) selected by cross-correlating 
the Gil de Paz et al. (2003) sample with the Sloan Digital Sky Survey Data Release 6 
(SDSS DR6). The sample includes 31 BCDs, which span a large range in galactic pa- 
rameters. Using a stellar population synthesis method, we derive the stellar populations 
and reconstruct the star formation histories for these BCDs. Our studies confirm that 
BCDs are not young systems experiencing their first star formation but old systems 
undergoing a starburst activity. The stellar mass-weighted ages can be as old as 10 Gyr 
while the luminosity-weighted ages might be up to ~ 3 orders of magnitude younger 
(~ 10 Myr) for most galaxies. 

Based on multi-wavelength data, we also study the integrated star formation prop- 
erties. The SFR for our sample galaxies spans nearly 3 orders of magnitude, from a few 
10~ 3 to ~ IMq yr" 1 , with the median value of ~ 0.1M Q yr _1 . We find that about 
90% BCDs in our sample have their birthrate parameter (the ratio of the current SFR 
to the averaged past SFR) b > 2 — 3. We further discuss correlations of the current 
SFR with the integrated galactic stellar mass and explore the connection between SFR 
and metallicity. 

Subject headings: galaxies: dwarf - galaxies: starburst - galaxies: stellar content 

1. Introduction 

Dwarf galaxies play an important role in our understanding of the formation and evolution of 
galaxies. In the hierarchical model of galaxy formation, they are proposed to be the building-blocks 
from which larger systems have been created by merging (Kauffmann et al. 1993). However, these 
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building-block galaxies are too small and faint to be studied at high redshifts. Therefore, of special 
interest are the studies aimed at the search for some local example galaxies which might have much 
in common with the galaxies at high redshifts. 

Blue compact dwarf galaxies (BCDs) are such local example galaxies, whose characteristics are 
believed to have been common among unevolved low-mass galaxies at intermediate to high redshift. 
BCDs are small, gas-rich (Hi mass fraction typically higher than 30%; e.g. Thuan & Martin 1981; 
Salzer et al. 2002) and metal-poor (1/50 < Z < 1/3 Z@; Hunter & Hoffman 1999) extragalactic 
objects. They have dramatically different properties compared to normal dwarf galaxies (Zwicky 
1966; Gil de Paz et al. 2003, hereafter G03), and are spectroscopically characterized by a faint, 
blue optical continuum accompanied, in most cases, by strong narrow emission lines, due to the 
intense star formation activity in one or several star- forming regions (Cairos et al. 2001). Stars are 
formed at high rates in BCDs (Fanelli et al. 1988), exhausting their gas content with a timescale 
of 10 9 yr (e.g. Lee et al. 2002), which is much shorter than the age of the Universe. This fact, 
combined with the low metal abundances, led Sargent & Searle (1970) to suggest that BCDs are 
primeval galaxies undergoing star formation for the first time. 

The question as to whether BCDs are truly young galaxies or they are old galaxies exhibiting 
a strong starburst has been widely discussed for many years (e.g. Sargent & Searle 1970; Searle et 
al. 1973; Schulte-Ladbeck et al. 2000, 2001a; Izotov & Thuan 1999, 2004; Aloisi et al. 2005). Most 
BCDs are known to have a red low surface brightness background of presumably older stars (e.g., 
Loose &: Thuan 1986; Gil de Paz &; Madore 2005; Caon et al. 2005). Moreover, some individual 
red giant stars (RGBs) in the nearest BCDs can be resolved (e.g., Lynds et al. 1998; Drozdovsky 
et al. 2001; Corbin et al. 2008). This evidence ruled out the possibility that the majority of BCDs 
began to form stars within the last billion years. 

An alternative view proposed by Searle et al. (1973) is that BCDs are chemically primitive 
objects which experience an episodic star formation history (SFH). The star formation in BCDs 
occurs in intense bursts which are interleaved by long quiescent periods. This traditional picture also 
seems to be problematic. Several nearby BCDs with resolved stars do not show evidence of long gaps 
(> 1 Gyr) in their recent star formation history (Schulte-Ladbeck et al. 2001a; Crone et al. 2002; 
Annibali et al. 2003; McQuinn et al. 2009). Sanchez Almeida et al. (2008) identified a large sample 
of quiescent counterparts to BCDs (QBCD) in SDSS Data Release 6 (DR6; Adelman-McCarthy 
et al. 2008), which might support this recursive BCD phase scenario. They also argued that the 
quiescent phase could last 30 times longer than the starburst phase, and thus the quiescent period 
of BCDs might be several hundred Myrs under the assumption of a single 10 Myr long starburst 
per BCD phase. However, Sanchez Almeida et al. (2008) showed that only ~ 15% of their QBCD 
candidates exhibit Ha absorptions and therefore are lack of star- forming activities. Combining 
the results described above with the fact that very few late-type and Hi selected galaxies have no 
observable Ha emissions (Meurer et al. 2006; James et al. 2008; Kennicutt et al. 2008), the SFR 
of the quiescent phase is probably not zero, but rather may be characterized by a continuous, low 
level of activity as discussed in Lee et al. (2009). 
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The detailed studies of SFHs of BCDs are rather difficult. For the most nearby systems, where 
individual stars can be resolved, the Color-Magnitude Diagrams (CMDs) synthesis method could 
be used to reconstruct the SFHs (e.g. I Zw 18, Aloisi et al. 1999, 2007; Mrk 178, Schulte-Ladbeck 
et al. 2000; I Zw 36, Schulte-Ladbeck et al. 2001a; NGC 1705, Annibali et al. 2003). However, 
this method can not be applied to more distant objects, which are the great majority of BCDs. 
The only way to study their stellar populations is to compare their integrated properties with the 
predictions of evolutionary synthesis models. There have been several works using (evolutionary) 
synthesis models to study the observed spectral energy distribution of BCDs (e.g. Mas-Hess & 
Kunth 1999; Guseva et al. 2001; Corbin et al. 2006). However, these studies only include one 
or a few galaxies. In this paper, we present a detailed study on the stellar populations of a BCD 
sample containing 31 galaxies, using a Simple Stellar Population (SSP) synthesis method. This is 
capable of yielding the various stellar components, the internal reddening (both stellar continuum 
and nebular line emission) and the SFH. The pure emission-line spectra allow us to measure the 
fluxes of emission lines more accurately. We can also study the star formation properties of these 
31 BCDs by the using the continuum-subtracted, narrow band Ha imaging data (G03) combined 
with the results presented in this work. 

The paper is organized as follows: Section 2 describes the BCD sample and our data reductions. 
Our results and analysis are given in Section 3. We discuss the uncertainties of the stellar population 
synthesis and the aperture effects of the fiber spectra in Section 4, and summarize our results in 
the last section. Where required we adopt a Hubble constant of Hq = 70 kms" 1 Mpc^ 1 , 17m = 0.3 
and f^A = 0.7. 

2. Sample, Data and Method 

2.1. Sample and Data 

There have been a number of criteria for selecting BCDs in literature. These criteria are 
commonly based on the galaxy's morphological properties and its luminosity (Zwicky & Zwicky 
1971; Thuan & Martin 1981), although some definitions are based on their spectroscopic properties 
(Gallego et al. 1996). BCDs classified by these different criteria may be confused with each other, 
e.g. a spectroscopically classified BCD may not be a BCD based on its morphological or luminosity 
properties (see G03). 

Using a unified concept of BCD galaxy, G03 compiled a large but not complete sample of BCDs 
(including 105 members with heliocentric recession velocities less than 4000 km s _1 ) from several 
exploratory studies (Haro 1956; Zwicky & Zwicky 1971; Smith et al. 1976; Kunth et al. 1981; 
MacAlpine & Williams 1981; Markarian et al. 1981, 1986; Zamorano et al. 1994, 1996; Alonso 
et al. 1999; Ugryumov et al. 1999), in which the galaxies were selected or discovered by (1) the 
ultraviolet excess, or (2) the presence of emission lines, or (3) the compactness and blue colors. The 
galaxies in the G03 sample were selected by putting forward a new set of quantitative classification 
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criteria, namely a combination of the galaxy's color, morphology and luminosity. As presented in 
G03, these criteria are 

/^_B,peak — ^i?,pcak < 1, 

MB,peak < 22 mag arcsec -2 , 
M K > -21 mag. 

where ^B iPea k and /x_R )Pea k are the peak surface brightness in B- and i?-band, respectively. 

By cross-correlating the BCD sample in G03 with the SDSS DR6, we obtained the optical 
spectra for 31 BCDs. In Table 1 we list some basic parameters for these 31 galaxies. The distances 
are all adopted from G03. Except for three galaxies whose distances have been determined by 
measuring the magnitude of the tip of the RGB (Mrk 178, Schulte-Ladbeck et al. 2000) and by 
using the period-luminosity relation of Cepheids (Haro 8 and VCC 0130, Macri et al. 1999), the 
rest were computed using the galactic standard of rest velocity of the galaxies. The gas phase 
oxygen abundances, both derived with direct and empirical methods, are taken from Zhao et al. 
(2010; hereafter ZGG10; see also the references therein). The direct method is to use electron 
temperature-sensitive lines (the so-called T e -method), whereas the empirical method (e.g. N2- 
method; Denicolo et al. 2002) measures the oxygen abundance through empirical relation between 
emission line ratio(s) and metallicity. As pointed out in Denicolo et al. (2002) and discussed in 
ZGG10 (also see Pettini & Pagel 2004 and Salzer et al. 2005), the metallicity measured with the 
N2-method is less precise than with the T e -method, and has an accuracy of ~ 0.2 dex. 

Therefore, our sample is not complete in any sense, which is common in all BCD studies in the 
literature, as it is difficult to achieve a complete sample for faint galaxies (however, see recent efforts 
to construct complete sample of dwarf galaxies from Blanton et al. (2005), Kennicutt et al. (2008) 
and Lee et al. (2009)). Fortunately, our sample is derived from the largest BCD sample to date and 
spans a large range in galactic parameters including integrated luminosity (Mb of —12.0 to —18.5), 
color (B — R of 0.0 — 1.5), central surface brightness and size (hr, o of 18.5 — 23.5 mag arcsec -2 and 
effective radius at i?-band of 0.09 — 2.2 kpc, see Gil de Paz & Madore 2005) of the underlying stellar 
population (USP), current star formation activity (0.03 — 1.6 Mq yr -1 kpc -2 ), oxygen abundance 
(12 + log(0/H) of ~ 1/20 - 1/2 Z Q ), and relative gas content (0.1 - 3.3 M q /L b ,q, for 24 galaxies 
that have Hi observations in the literature, see Zhao et al. (2010, in prep, and references therein). 
Because of the complexity and diversity of galaxies, we need a large sample to probe the full range 
in galactic parameters in the future work. 

In order to check whether this subsample is representative of the total sample of G03 with 
respect to the integrated (B — R) color and the absolute B magnitude (Mb), we plot the (B — R) 
and Mb distributions for the G03 sample and our sample using blank and hatched histograms, 
respectively, in Figure 1. We can see that our sample has similar distributions of (B — R) and Mb 
to the G03 sample, except that our sample contains a smaller fraction of galaxies with (B — R) 
bluer than 0.5. The B- and .R-band photometric data are directly taken from G03, where the 
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Table 1. Basic properties of the BCD sample 



Distance 12 + log (O/H) 3 







T 


N2 


HS 0822+3542 


10.1 


7.42 + 0.03 


7.47 


HS 1400+3927 


21.1 


8.11 + 0.06 


8.08 


HS 1440+4302 


38.0 


8.09 ± 0.03 


8.08 


HS 1609+4827 


42.6 


8.14 + 0.14 


8.37 


Haro 2 


21.6 


8.38 + 0.03 


8.64 


Haro 3 


14.4 


8.46 + 0.10 


8.39 


Haro 8 


16.0 


8.35 + 0.04 


8.16 


I Zw 123 


11.7 


8.07 + 0.02 


7.99 


II Zw 71 


18.7 


8.24 + 0.17 


8.42 


Mrk 1313 


31.8 


8.22 + 0.10 


8.05 


Mrk 1416 


33.8 


7.84 + 0.02 


7.94 


Mrk 1418 


11.4 




8.48 


Mrk 1423 


20.4 




8.52 


Mrk 1450 


14.7 


7.96 + 0.02 


7.93 


Mrk 1480 


27.2 


8.04 + 0.05 


8.03 


Mrk 1481 


27.3 




8.11 


Mrk 178 


4.2 


7.92 + 0.02 


7.88 


Mrk 409 


21.3 




8.80 


Mrk 475 


9.7 


7.93 + 0.02 


7.89 


Mrk 67 


14.3 


8.08 ± 0.08 


7.90 


SBS 1054+504 


20.1 


8.26 + 0.10 


8.34 


SBS 1147+520 


18.9 




8.19 


SBS 1428+457 


35.3 


8.42 + 0.05 


8.18 


SBS 1533+574 


49.6 


8.07 + 0.07 


8.11 


UGCA 184 .. . 


23.0 


8.03 + 0.01 


7.89 


UM 323 


28.3 


7.96 ± 0.04 


8.00 


UM 439 


14.0 


8.08 ± 0.03 


7.93 


UM 452 


19.5 


8.27 + 0.20 


8.40 


UM 456A .... 


24.5 




8.21 


UM 491 


27.6 




8.31 


VCC 0130 .... 


16.0 




8.28 



a Oxygcn abundances derived with the T e -mcthod 
and N2-method, respectively, see ZGG10 for more 
details. 
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foreground Galactic extinctions have been corrected using Ab values determined following Schlegel 
et al. (1998) and the Galactic extinction law of Cardelli et al. (1989). 
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Fig. 1. — Number distribution of the BCDs in (B — R) color (upper panel) and absolute B 
magnitude (bottom panel). Blank histogram: G03 sample; Hatched histogram: our sample. 

The purpose of this work is to constrain the SFHs of these BCDs. However, as shown in Figure 
2, the BCDs in our sample are all extended objects. The SDSS spectra are from the starlight that 
found its way down a 3"-diameter fiber (with a physical size of ~ 60 pc at the nearest object Mrk 
178 and ~ 720 pc at the most distant object SBS 1533+574) that have been positioned on bright 
Hn regions for these BCDs (see Figure 2 for a direct view). Therefore, we need to check to what 
extent the SDSS fiber spectrum represents the light of the whole galaxy. 

In Table 2, we list two quantities related to the covering fraction of the fixed-size aperture 
for all BCDs. One is the ratio (r^^/petroRgo) between the fiber radius (l"-5) and the Pertrosian 
radius containing 90 percent of the galaxy light in r-band, and the other is the light fraction, f r , 
captured by the fiber. The f r is calculated using the total galaxy magnitude and the magnitude 
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Fig. 2. — SDSS r-band images of the BCDs. The overlaid circle in each panel shows the position 
of the SDSS 3"-diameter fiber. The sizes of all images are 40" x 40". 
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Fig. 2. — Continued 

inside the fiber, explicitly f r = 100 x l0 _a4 (' mflbc '' _ ' mpcrtr °). The median values of rgbcr/petr 0-^90 
and f r are 0.37 and 14.4%, respectively. For our sample galaxies, therefore, the SDSS spectrum 
could have much bluer (g — r) color and less prominent 4000A break (D n (4000)) (see below) than 
the integral spectrum. 

In Figure 3, we plot the observed SDSS spectra (black line) together with the integrated ugriz 
broad-band spectral energy distributions (SEDs; red dashed line) for all galaxies in our sample. 
The flux has been normalized to the i-band value. For most galaxies, the continuum shape of the 
fiber spectrum longward of r-band is very consistent with that of the broad-band SED, whereas 
the spectrum shortward of r-band starts to deviate from the broad-band SED to varying degrees. 
It suggests that the aperture effects might much affect our fitting results presented in the next 
section for galaxies whose fiber spectra apparently differ from the broad-band SEDs blueward of 
r-band. Fortunately, ten out of the 31 galaxies have been observed using drift-scanning technique 
(Moustakas & Kennicutt 2006, hereafter MK06; private communication). The integrated spectra 
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for these galaxies are superposed in Figure 3 with blue lines. Again, the aperture effects, especially 
of the 4000 A break, are clearly demonstrated. We will return to the question of how much the 
aperture effects will affect our results in the next and the discussion sections. 

Kewley et al. (2005) have shown that the fiber needs to capture more than 20 percent of 
the galaxy light to minimize the aperture effects in the spectral measurements. According to this 
threshold we divide our galaxies into three subsamples: f r < 20%, 20% ^ f r < 100% and f r ~ 100% 
(integrated spectra), and in the following we refer them to SI, S2 and S3, respectively. There are 
16, 5 and 10 galaxies in SI, S2 and S3, respectively. 



2.2. Stellar Population Synthesis 

In order to obtain the stellar populations for these BCDs and to correct for the stellar ab- 
sorption in the emission lines, we here model the stellar contribution in the SDSS spectra through 
the modified version of the stellar population synthesis code, STARLIGHT0 (Cid Fernandes et al. 
2005 (C05), 2007; Mateus et al. 2006; Asari et al. 2007), which adopted the stellar library from 
Bruzual 8z Chariot (2003). The code does a search for the linear combination of SSPs to match a 
given observed spectrum (0\). The model spectrum (M\) is: 

M x = M\ J2f=i Xjb j:X rx ® 

where bj^x = L^ SP (tj, Zj)/Lf^ p (tj, Zj) is the spectrum of the j th SSP normalized at Ao, rx = 
-^q-oa(A x -Ax ) j g reddening term, x is the population vector, Mx is the synthetic flux at the 
normalization wavelength, A* is the total number of SSPs, G(t>*, cr*) is the line-of-sight stellar 
velocity distribution, modeled as a Gaussian centered at velocity and broadened by a*. The 
match between model and observed spectra is calculated by 

X 2 (x, Mx ,A v ,v+, a,) = J2x=i 1(0 x ~ M x ) w x f, 

where the weight spectrum wx is defined as the inverse of the noise in Ox- For more details, 
please refer to C05 and Mateus et al. (2006). The SSP library follows the work of SEAGal 
Group, and is made up of N* = 100, including 25 ages (from 1 Myr to 18 Gyr) and 4 metallicities 
[Z = 0.005, 0.02, 0.2, and OAZq, for our metal-poor sample). The spectra were computed with 
the Salpeter (1955) initial mass function (IMF), Padova 1994 models and the STELIB library (Le 
Borgne et al. 2003). The intrinsic reddening is modeled by the foreground dust model, using the 
extinction law of Calzetti et al. (1994) with Ry = 4.05 (Calzetti et al. 2000) for these star-forming 
galaxies. 



STARLIGHT & SEAGal: 



http: / /www. starlight, ufsabr/] 
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Fig. 3. — The observed SDSS spectra (black line) and broad-band SEDs (ugriz; red dashed line) 
for these 31 BCDs. We also overplot the integrated spectra from MK06 (blue line) for 10 BCDs. 
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Fig. 4. — Results of the spectral fitting for Haro 2. The top left panel shows the logarithm of 
the observed (F°; black line) and the synthetic {F^ 1 '; red line) spectra. The F® — F^ 1 residual 
spectrum is shown in the bottom left panel. Spectral regions actually used in the synthesis are 
plotted with a black line, while masked regions are plotted with a green line. Panels in the right 
show the population vector binned in the 25 ages of SSPs used in the model library. The top right 
panel corresponds to the population vector in flux fraction, normalized to Ao = 4020 A, while the 
corresponding mass fractions vector is shown in the bottom right panel. 



The SDSS spectra cover 3800-9200 A, with a resolution (A/AA) of 1800 < R < 2100 and 
sampling of 2.4 pixels per resolution element. The fiber used in the SDSS spectroscopic observations 
has a diameter of 3" on the sky. Prior to the synthesis, the Galactic extinction has been corrected 
with a combination of the extinction law of Cardelli et al. (1989) and the Ab value from Schlegel 
et al. (1998) as listed in The spectra are transformed into the rest frame using the redshifts 

given in the FITS header. The SSPs are normalized at Ao = 4020 A, while the observed spectra 
are normalized to the median flux between 4010 and 4060 A. The signal-to-noise ratio (S/N) is 
measured in the relatively clean window between 4730 and 4780 A. Masks of 20 — 30 A around 
obvious emission lines are constructed for each object individually, and more weights are given to 
the strongest stellar absorption features such as Call K A 3934, and the Can triplets, that are less 
affected by nearby emission lines. 



2 http://ned www.ipac.caltech.edu/ 
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For our sample, the S/N varies between 6.8 and 72.0 (see Table 2), and most of them are 
larger than 20. Generally, the fitting results for high S/N objects are better than those for low S/N 
ones. Inspecting the fitting results, we find that the goodness of fitting (x 2 value) also somewhat 
depends on the absorption line equivalent widths (e.g. EW of Can K). A typical example of our 
fitting result for Haro 2 is shown in Figure 4. After subtracting the best-fit model spectrum from 
the observed one, we obtain the pure emission-line spectrum as shown in the bottom left panel of 
Figure 4, where we can derive accurate fluxes for all the emission lines with the onedspec.splot task 
in IRAP0. 

3. Results and analysis 

3.1. Mean Stellar Age 

In the top panels of Figure 5, we show the distributions of the mean stellar ages (C05) estimated 
for all of the three subsamples (Figure 5a, 56 and 5c for SI, S2 and S3, respectively). They are the 
mass-weighted mean stellar age, 

(log t*) M = J2vj lo S *i» (!) 

and the light-weighted mean stellar age, 

(log t+) L = ^2,Xj log tj. (2) 

where fj,j and Xj represent the fractional contributions to the stellar mass and luminosity of the 
SSP with age tj respectively. N* is the number of SSPs. Based on these two definitions, it is 
easy to understand that (t*) M represents the epoch of the stellar population formation and these 
stars nowadays contribute most significantly to the stellar mass in galaxies. Consequently, {t*) M is 
associated with the mass assembly history of galaxies. On the other hand, (£*)l renec ts the epoch 
of the formation of massive and bright O and B stars, frequently associated with recent starburst 
activities. Therefore, (t*)^ is strongly affected by the recent SFH of a given galaxy. The fitting 
results of these two kinds of ages for our BCD sample are summarized in Table 2. According to 
C05 and Mateus et al. (2006), the uncertainties of these two parameters depend on the S/N of the 
input spectra. In general, the rms of the fitted (log t*) M is ~ 0.2 dex for S/N < 10 and ~ 0.1 dex 
for S/N > 10, while the rms of the fitted (log t*) L is ~ 0.15 dex for S/N ~ 5 and < 0.1 dex for 
S/N > 10. 



3 IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under a cooperative agreement with the National Science 
Foundation. 
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Table 2. 


Derived and fitting 


results 


for the BCD sample 




petroi?go b 




fr C 




(1°S Ml 


< lo g **>M 


Galaxy 


(arcsec) 


-Kfibcr/-R9l 


, (%) 


S/N d 


(yr) 


(yr) 


HS 0822+3542 a 


1 9 


0.80 


49.5 


6 8 


7.68 


8.09 










2 6 


7.10 


9.71 


HS 1400+3927 . 


3 


0.51 


25.3 


22.7 


7.40 


8.27 


HS 1440+4302 . 


2 9 


0.52 


26.7 


12.2 


7.43 


9.85 


HS 1609+4827 . 


3 9 


0.39 


13.6 


34.4 


7.23 


9.87 


Haro 2 a 


5 8 


0.26 


13.4 


72.1 


7.36 


9.94 










62.2 


7.76 


9.91 


Haro 3 a 


7.4 


0.20 


5.0 


52.2 


7.76 


9.56 










31.0 


7.99 


9.94 


Haro 8 a 


8 9 


0.17 


7.2 


25.9 


7.47 


10.02 










21.8 


8.32 


9.89 


I Zw 123 a 


3.7 


0.41 


17.2 


12.1 


6.99 


9.98 










34.8 


7.19 


9.53 


II Zw 71 


S 8 


0.17 


1.3 


46.0 


7.89 


9.34 


Mrk 1313 


3 1 


0.48 


19.3 


36.4 


7.68 


8.92 


Mrk 1416 


4 5 


0.34 


15.3 


25.7 


6.87 


7.10 


Mrk 1418 a 


7.0 


0.21 


3.7 


49.0 


7.96 


8.86 










49.6 


7.94 


9.90 


Mrk 1423 


11.2 


0.13 


6.4 


33.5 


7.75 


8.67 


Mrk 1450 a 


2.1 


0.70 


33.6 


20.5 


7.08 


9.96 










17.7 


7.03 


9.74 


Mrk 1480 


3.2 


0.47 


22.6 


31.4 


7.01 


9.58 


Mrk 1481 


6 6 


0.23 


5.8 


7 7 


7.37 


9.87 


Mrk 178 a 


8 4 


0.18 


7.3 


34.0 
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a The second line gives the fitting results based on the integrated spectra. 
b The Petrosian radius containing 90 percent of the galaxy flux in r-band. 
c The fraction of the r-band flux covered by the SDSS fiber. 
d The signal-to-noise in the S/N window. 
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In each panel of the top row of Figure 5, the solid line shows (log£*} A j, while the dashed line 
shows (log i*)r- For all of the three subsamples, (logt*)^ is in the range of several Myr to a few 
100 Myr, with a trend of the median (logi*)^ towards larger value from SI to S3. 




a (log U) » (log t.) a (log t.) 



Fig. 5. — Statistics on the fitting results (from left to right for subsamples: SI, S2, and S3). In 
each panel, the mass-weighted and light-weighted results are shown with solid and dashed lines, 
respectively. The filled and (thick) blank histograms in the right panels (c and /) show the results 
derived with the fiber and integrated spectra, respectively. Top panels: Number distributions of the 
BCDs in mean stellar ages. Bottom panels: Number distributions of the BCDs in the dispersions 
of (log**)- 

As shown in Figure 5a, for the SI sample, whose spectra might severely under-sample the 
member galaxies, (logi*)jy ranges from ~ 10 Myr to ~ 10 Gyr, and 11 out 16 galaxies have 
(logt*) A/ > 1 Gyr. For the S2 sample, whose spectra may be representative of the member galaxies, 
(logi*) A/ is in the range of several 100 Myr to 10 Gyr (Figure 56). For the S3 sample, which have 
integrated spectra, all of (logi*) M are much larger than 1 Gyr (blank, thick histogram in Figure 
5c). These results indicate that the aperture bias will lead to an underestimate of (logi*)^, 
which can also be seen from Figure 5c, in which the filled histogram shows the results derived 
from the fiber spectra. To check the aperture effect further, we compare results derived from the 
integrated spectra with those from the fiber spectra for the 4 galaxies having (logi*) A/ < 1 Gyr, 
and find that the mean values of A (logi*}^ (= {^ogt^) Lintc ^ iatcd — { l °SU) L 3 ») and A (logt*) M (= 
(logt+) M intcgratcd — {logt+) M 3 //) are 0.34 ± 0.49 and 0.90 ± 1.38, respectively. Thus, while both 
(logi*) L and (logt*) A/ are underestimated based on the fiber spectra, (logt*) A/ suffers much more 
severely from the aperture effects. 
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In our sample, Mrk 178 is the only galaxy which has been studied extensively using CMDs in 
literature. Schulte-Ladbeck et al. (2000) found RGB stars with age of at least 1 Gyr in this galaxy 
and concluded that this galaxy should be an old galaxy. The (logt*) M derived from the integrated 
spectrum is as old as ~ 5 Gyr, and is well consistent with the results of Schulte-Ladbeck et al. 
(2000). 

In summary, although a part of our sample may suffer from the aperture effects, our results 
are generally consistent with the findings through deep photometric observations for nearby BCDs 
(e.g. Schulte-Ladbeck et al. 2001a; Aloisi et al. 2007), and clearly show that BCDs are not truly 
primordial galaxies, rather, they are old systems undergoing transient periods of starburst. 

To investigate the SFH of BCDs in more detail, we calculated the light-weighted and mass- 
weighted standard deviations of the log age, which are defined as the following (C05), 

1/2 



0L(log U 



IV* 



^2 X j (log tj - (logt 



*/L) 



(3) 



and 



1/2 



a M (1ogU)= ^fijQogtj - (logi*) M ) 2 . (4) 
J= l 

These two higher moments of the age distribution could be used to distinguish galaxies dominated 
by a single population from those which had bursty or continuous SFHs. 

The bottom panels in Figure 5 displays the number frequency histograms of BCDs in these 
two parameters. Both of the median values of cr^f (log t±) and ox (log £*) are significantly larger than 
zero, indicating that BCDs are not dominated by a single population. ox(log t+) has a larger median 
value than <7jy/(log t*). This suggests that the majority of the stellar mass in BCDs is assembled in a 
relatively narrow time range, while the star formation activities might occur repeatedly during the 
lifetime of a BCD. At the same time, we can see from Figure 5/ that the aperture effects generally 
result in an underestimate of ox(logi*) and an overestimate of <7M(logt*). 



3.2. Luminosity and mass fractions of different populations 

A coarse but robust description of the SFH of a galaxy may be obtained by binning x (ju) into 
'young' (tj < 10 8 yr), 'intermediate-age' (10 8 < tj < 10 9 yr), and 'old' (tj > 10 9 yr) components (xy 
(/iy), x\ (m), and xo (a*o)j respectively). As shown in Figure 6, our stellar synthesis result indicates 
that the luminosity and stellar mass of the BCDs are contributed by different components. For 
most galaxies, the young population contributes the largest fraction of the 4020 A monochromatic 
radiation, while the old population contributes most stellar mass. 

For the S3 sample, more than 50% of the 4020 A monochromatic radiation comes from the 
young population for 8 out of 10 galaxies, while more than 80% of the stellar mass is contributed 
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Fig. 6. — Number distributions of the BCDs in luminosity (top half page) and mass (bottom half 
page) fractions for the young (t < 10 8 yr), intermediate (10 8 < t < 10 9 yr) and old (t > 10 9 
yr) stellar populations. Solid and dashed histograms show the integrated and fiber spectra-based 
results, respectively. For each of the half page, from top to bottom: SI, S2, and S3. 
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by the old population for all of the galaxies. For the SI and S2 samples, which might suffer severely 
from aperture effects, the young and inter-mediate populations have more contributions to the 
luminosity and mass fractions than the S3 sample. About 1/3 — 1/2 galaxies in SI and S3 have more 
than 20% of their stellar mass and light contributed by the intermediate-age population. Whereas 
for the S3 sample, the intermediate-age population only contributes < 20% light and stellar mass 
for all of the galaxies (only one exception that the intermediate-age population contributes 30-40% 
light). 

Our results (specially from the integrated spectra) are in agreement with the scenario that 
BCDs are old galaxies with episodic starbursting activities. The scarcity of intermediate-age com- 
ponent suggests that BCDs should not have a SFH with a constant (current) SFR. However, we 
can not discriminate the scenario that there exists extensively long quiescent phases during a BCD 
galaxy lifetime from that a BCD is still forming stars with a very low SFR during its "quiescent" 
periods. 



3.3. Extinction and the Age-Extinction Degeneracy 

There are several independent ways to measure the internal extinctions. In this paper, we 
use the continuum and nebular lines to make separate estimates of the internal extinctions in the 
BCDs. starlight also returns an estimate of the stellar visual extinction, Ay, modeled as due to 
a foreground dust screen. 

By making the reasonable assumption that the intrinsic Balmer line ratios are equal to the 
Case B recombination and using a Calzetti et al. (1994) reddening law, we could estimate the 
nebular extinctions from the observed Ha/H/3 and H7/H/3 Balmer decrements, which are 

A v a = 7.93 x log(^M) (5) 



Ay 1 = -19.04 x Jog(^i^) (6) 



lH7_ mn^ ,„„/V%- 

In*// lap 

where Fh q / F^p ,Fh 7 / F^p and Ina/Inp, In^/Iup are the observed and intrinsic Balmer decrements, 
respectively. Here we adopt the intrinsic ratios of Iua/Inp and In^/Iup to be 2.86 and 0.469 
(Brocklehurst 1971), respectively, for an electron temperature of 10 4 K and an electron density of 
100 cm" 3 . We prefer using equation (5) to calculate the extinction since H7 has relatively lower S/N 
than Ha and equation (6) is much more sensitive to the measurement errors compared to equation 
(5). The observed fluxes of Ha and H/3 are measured through the "pure-emission" spectra, and 
are presented in Table 3. The errors in the line fluxes are calculated according to files attached to 
each spectrum which contain the error in each pixel. Additionally, an error of ~ 2% in the emission 
line fluxes should be accounted for the uncertainties of standard star absolute fluxes used for flux 
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calibration (e.g., Oke 1990). For sources with saturated Ha, the adopted extinction is calculated 
using equation (6). The Ay values calculated through Balmer decrements are given in Table 3. 

We compare these two extinctions for our BCD sample, and find that Ay, which is returned 
by starlight, in general, is comparable with Ay^neb- The median value of Ay^^/Ay is about 1.0, 
which is different from that found in the detailed studies of nearby star forming galaxies by Calzetti 
et al. (1994), who postulate that nebular line emission is attenuated by roughly twice as much dust 
as the stellar continuum (see also Chariot & Fall 2000). This result may be due to the well known 
age-extinction degeneracy (e.g. Gordon et al. 1997), which acts in the sense of confusing old, less 
dusty systems with young, dusty ones and vice versa. 

The age-extinction degeneracy will affect the following fitted parameters, e.g. the luminosity- 
and mass-weighted mean stellar ages, the luminosity and mass fractions of different populations, 
and the stellar mass, to varying degrees, and we will further discuss the uncertainties caused by 
this effect in Section 4.1. 



Through the spectral synthesis, starlight presents the current stellar mass, M*, which corre- 
spond to the light entering the fiber. The total stellar masses of the galaxies are obtained from M* 
divided by /, where / is the fraction of the total galaxy luminosity in the i-band inside the fiber. 
For most galaxies, this can give a good estimate of the total stellar mass. However, this method 
will significantly underestimate the stellar mass for galaxies whose fiber spectra badly represent the 
whole galaxy in our sample. Therefore, for galaxies having integrated spectra we prefer to using the 
stellar masses derived from the integrated spectra. For galaxies whose Ay neb are significantly over- 
estimated and which do not have integrated spectra, the stellar masses are derived by re-running 
starlight with the same parameters as used in Section 2.2 except that we limit Ay to the value 
of Ay neb ( see Section 4.1). After these procedures, the stellar mass of three galaxies (Mrk 1416, 
Mrk 1423 and SBS1147+520) might be still much underestimated since their fitted (logi*)^ are all 
a bit less than 1 Gyr. We remove these three galaxies from further analysis involving stellar mass. 



Using the relation between SFR and nebular emission line (Ha) luminosity (Kennicutt 1998), 
namely 



3.4. 



Stellar Mass 



3.5. Star formation properties 



SFR(M year" 1 ) = 7.9 x 10" 42 L(Ha) (ergss -1 ) 



(7) 
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Table 3. Observed and derived properties for the BCD sample 
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a In units of 10 -17 erg s _1 cm~ 2 . 

b Internal extinction derived with the flux ratio of Balmer lines. 

c Internal extinction-corrected luminosity of Ha, in units of 10 38 erg s _1 . For spectra that the Ha lines are saturated, 
we use 2.86xL H o to estimate the Ha luminosity. 

d Star formation rate, derived with the Ha luminosity in Col (5) 

e The ratio (expressed as a percentage) of our Ha luminosity to G03's 

f Total star formation rate, calculated based on the internal extinction-corrected Ha luminosity given in G03. 
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Fig. 7. — Panel (a) shows the correlation of SFR surface density between £sfr 3// and £sfr s - The 
other three panels plot the number distributions of the BCDs in (b) fractions of the Ha (solid) 
and r-band (dashed) fluxes inside the fiber, (c) star formation rate, and (d) birthrate parameter, 
b. The solid line in panel (a) is a reference line for the case when the two quantities are the same. 
The arrow in each panel marks the position of the median value. 

we can estimate SFR for these BCDs. Prior to the conversion, the internal extinction is corrected by 
using the attenuations calculated in the previous section using the Balmer line ratios and assuming 
the Calzetti et al. (1994) reddening law. The Ha luminosity and SFR covered by the SDSS fiber 
are listed in Table 3, Columns 5 and 6, respectively. 

We also used the Ha luminosities measured by G03, after the correction of the internal redden- 
ing using the Ay, neb values measured here, to calculate the integrated SFRs of the BCD sample. 
This is because that the G03 measurements of the Ha luminosity are based upon narrowband 
imaging (having been corrected for the [N n] contribution; please see the Appendix in G03) and do 
not suffer from the aperture effects that plague the spectral fluxes. The integrated SFRs are given 
in Table 3, Column 8. In the following, parameters or statistics relating to SFR are all based on 
the integrated SFRs unless otherwise stated. 

In Figure 7 we demonstrate the integrated and 3"-sized star formation properties of these 31 
BCDs. Figure 7a illustrates the correlation between the SFR surface densities. The x-axis, SsFR e , 
is the integrated SFR normalized to the area defined by the scale length measured from exponential 
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fits to i?-band images (Gil de Paz & Madore 2005). The y-axis, £sfr 3 /m is the SFR surface density 
measured within the SDSS 3"-diameter aperture. On one hand, these two parameters show an 
expected correlation that a galaxy with higher Esfr 6 also has higher £sfr 3 " ■ On ^ ne other hand, 
Ssfr 3 „ is larger than SgFR e for most galaxies in this BCD sample, which is due to that these BCDs 
are extended objects while the SDSS fiber targets at the brightest part in most cases. As shown 
by the solid histogram in Figure 76, the fraction of the Ha flux (Rft a ; see Table 3) captured in 
the fiber spectrum is small for most galaxies, with a median value of 35.4%. For nearly one fourth 
of the galaxies, i?Ho is less than ~ 10%. However, the median value of Ru a is much larger than 
that of R r {= f r ) (see the arrows in Figure 76), which indicates that the distribution of Ha flux is 
relatively more concentrated. 

Figure 7c shows the number distribution of the SFRs for these 31 BCDs. More than half 
galaxies in this sample have SFRs less than 0.1 M® yr _1 . The median value of the SFRs is 0.093 
Mq yr _1 . Using the stellar mass derived in Section 3.4, we can determine the birthrate parameter 
(Kennicutt 1983), 6, for this sample of BCDs. The parameter 6 is defined as the ratio of the current 
SFR to the average past value, SFR/ < SFR > pas t, where <SFR> past is defined as the stellar mass, 
M*, divided by the stellar age, t. To minimize the aperture effect, the stellar age we adopted here 
is the age of the oldest population given by the fitting procedure. 

Figure 7c? displays the distribution of log 6. The 6 parameter is usually used to separate 
starburst galaxies from normal star forming galaxies. A galaxy with 6 > 2 ~ 3 is defined as a 
starburst galaxy (e.g., Hunter & Gallagher 1986; Salzer 1989; Gallagher 2005; Brinchmann et al. 
2004; Kennicutt et al. 2005). For the subsample, the median value of log 6 is about 0.9 and more 
than 90 percent of the BCDs (26/28) have log 6 > 0.3, which indicates an ongoing starburst activity. 

Schulte-Ladbeck et al. (2001b) claimed that there might exist a correlation between the oxygen 
abundance and SFR (see also Hopkins et al. 2002). However, as shown in Figures 8a and 86, we 
find no evidence of obvious correlation between these two parameters. A Spearman rank correlation 
coefficient analysis indicates a very weak correlation (p = 0.25 and 0.38 at a 1.4<7 and 1.8<7 level 
of significance for Figures 8a and 86, respectively). Therefore, there is little correlation, if any, 
between SFR and oxygen abundance, which is supported by the fact that there exist large scatters 
in the plots in Schulte-Ladbeck et al. (2001b; Figure 4) and Hopkins et al. (2002; Figure 6). This 
argument may also be supported by the results obtained by Kobulnicky & Skillman (1996, 1997), 
in which the authors show that their measurements of individual Hll regions in dwarf Irregular 
galaxies suggest localized chemical enrichment from the active star formation generally does not 
occur. 

We can gain further insights on the nature of star formation for BCDs by examining the 
correlation between galaxy stellar mass and SFR. Noeske et al. (2007), Elbaz et al. (2007), Daddi 
et al. (2007) and Pannella et al. (2009) have shown that star formation and stellar mass define a 
tight correlation in galaxies at z ~ 1 and z ~ 2, with rough proportionality (logarithmic slope of 
~ 1.0). Similar proportionality is also seen at z = in data from the SDSS (Elbaz et al. 2007). 
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Fig. 8. — Correlations between Ha-based SFR and (a) oxygen abundance derived with the N2- 
method, (b) oxygen abundance derived with the T e -method, and (c) stellar mass. There are obvious 
trends for these BCDs with higher stellar mass to also show higher SFRs. However, no clear 
correlation between metallicity and SFR can be found. 

We illustrate the relationship between SFR and stellar mass, for our BCDs in Figure 8c, which 
indicates a good correlation between SFR and M*. A Spearman rank correlation coefficient analysis 
gives p = 0.68 at a 3.5a level. A nonweighted least-squares linear fit, using a geometrical mean 
functional relationship (Isobe et al. 1990), to these 28 BCDs gives log (SFR) cx (1.21±0.18) log M*, 
which indicates that the SFR increases with stellar mass almost linearly. Our slope of the log(SFR)- 
log(M^) is in agreement with those found by Noeske et al. (2007), Elbaz et al. (2007), Daddi et 
al. (2007) and Pannella et al. (2009) for more massive local and high-redshift galaxies. However, 
there is a significant scatter in the trend, i.e., for a given stellar mass, the SFR may vary over an 
order of magnitude. 
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4. Discussion 

In this section we mainly discuss the uncertainties of the spectral synthesis and the aperture 
effects. 

4.1. Uncertainties of the Spectral Synthesis 

The SEAGal Group has used the starlight to analyze several large samples of SDSS galaxies 
as shown in their series of works (C05; Cid Fernandes et al. 2007; Mateus et al. 2006; Asari 
et al. 2007). As mentioned in Section 3.1, they have tested the uncertainties of the resulted 
stellar populations. In the study of C05, they found that the individual components of x are very 
uncertain, whereas the binned vectors of x, i.e. the young, inter-mediate and old populations, have 
uncertainties less than 0.05, 0.1 and 0.1, respectively, for S/N>10. 

For each fitting, starlight provides the last-chain-values of the contributed light (mass) 
fraction of each SSP, x 2 an d mass, for seven Markov chains. For the results of our sample, the 
median values of the rms of these adopted values are 1.8%, 17.4%, 8.8%, 8.1%, 21.0%, 1.9%, 
0.3%, 0.5%, 1.8%, 0.06%, 5.1% for x Y , x h x , /z Y , W, ^O, (logt*) L , (logU) M , Ay, X 2 and M*, 
respectively. Therefore, the uncertainties of the resulting stellar populations and stellar mass will 
not much affect our conclusions. 

Most of the age sensitivity comes from the 4000 A break and continuum shape in the fitting, 
thus degeneracy with extinction should be discussed carefully. Stellar mass-related parameters 
(e.g. (logt*) M and M*) will suffer more severely from this degeneracy because of the non-constant 
mass-to-light ratio of stars. As described in Section 3.3, the dust extinction is considered as a 
variable parameter in starlight, and it is overestimated (comparing to ^4y, n eb) for at least half of 
our sample, which might result in a large underestimation of (logi*)^-. For example, out of the 11 
galaxies whose fitted (log£*) M are less than 1 Gyr based on the fiber spectra (see Table 2), seven 
have their Ay overestimated. To check to what extent this overestimation will affect the fitted 
results, we re-run starlight using the same parameters except that we limit Ay to the value of 
Ay^neb f° r ten galaxies whose Ay /Ay, neb > 1-5. We find that A (logi*) L = 0.24±0.10, A (log^) M = 
1.15 ± 1.03, and AlogM* = 0.41 ± 0.27. Among these ten galaxies, six have (logt*}^ < 1 Gyr 
(as listed in Table 2), and we obtain A(logi*) L = 0.25 ± 0.12, A(logi*} M = 1.74 ± 0.92, and 
AlogM* = 0.56 ± 0.24. Therefore, the uncertainties of our fitting results are dominated by the 
age-extinction degeneracy. 

4.2. Aperture effects 

The SDSS is a fiber-based survey, and thus we should consider the biases introduced by the use 
of small apertures to measure the galaxy spectra. Kewley et al. (2005) have examined this problem 



in detail and they conclude that 20% of the galaxy light is required to minimize the aperture effects. 
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Fig. 9. — The mean light- (upper panels) and mass-weighted (bottom panels) stellar age as a 
function of the (u — g) and (g — r) colors within the 3"-diameter region. Open and solid circles 
show the spectra with f r < 20% and f r > 20%, respectively. 

As described earlier, the SDSS spectra for our sample galaxies only cover a small fraction of 
the galaxy light, and the median f r is 14.4%. The fiber spectra are bluer than the integrated 
SEDs for most galaxies. Based on the fiber spectra, in general, x\ is overestimated and xo is 
underestimated (see Figure 6), resulting in an underestimation of both (logt*) L and (logi*) M . 
Therefore we want to check to what extent our estimates of the stellar properties may be biased 
on the basis of the fiber spectra. In Figure 9 we plot the dependence of the mean stellar ages on 
the fiber (u — g) and (g — r) colors. We find that (logi*)^ correlates both with (u — g) and (g — r) 
(see Figures 9a and 96). Using these two correlations, we can make a coarse estimation of the 
bias in (logi*)^. For our sample, the 3"-regions are generally 0.15 and 0.17 mag (median value) 
bluer than the whole galaxies in (u — g) and (g — r), respectively, which result in 0.30 and 0.42 dex 
underestimations of (logi*) L respectively. However, no dependence of (logt*) A/ on these two colors 
can be established. This is because that {logt ir ) M has much less direct relation with the observed 
spectrum and is much more affected by the age-extinction degeneracy than (logi*)^ due to the 
non-constant mass-to-light ratio of stars. As shown in Section 3.1, the aperture effects generally 
cause an underestimation of (logt*) M for our sample galaxies, and thus, combining with the results 
derived with the integrated spectra, they might not much affect our conclusion that BCDs are old 
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galaxies experiencing starbursting activities. 

5. Summary 

In this paper we present the detailed results of our stellar population synthesis for an overlapped 
sample of 31 BCDs between G03 and SDSS DR6, using both integrated and fiber spectra. Our 
studies show that BCDs are not young systems experiencing their first star formation activities 
but old systems with episodic star formation histories, which confirm previous studies. For most 
galaxies, the stellar mass-weighted ages for the 3"-regions in BCDs are as old as 10 Gyr while the 
luminosity-weighted ages could be ~2 to 3 orders of magnitude younger. While for the 10 galaxies 
having integrated spectra, the stellar mass-weighted ages are all larger than several Gyrs. Old 
population contributes the majority of stellar mass while young population contributes most of 
light in the optical. 

The SFR of more than half sample galaxies is less than 0.1 Mq yr _1 , and the median SFR 
is 0.093 Mq yr^ 1 . For 28 BCDs which have their stellar mass measured reliably, we calculate the 
birthrate parameter, b, and find that about 90% galaxies have b > 2 — 3, which indicates that 
these galaxies are undergoing a starburst. The current star formation rate correlates well with the 
integrated galactic stellar mass. The nebular metallicity shows a very weak correlation with the 
current SFR. 
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